Helical Luttinger liquid in topological insulator nanowires 
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We derive and analyze the effective low-energy theory for interacting electrons in a cylindrical 
nanowire made of a strong topological insulator. Three different approaches provide a consistent 
picture for the band structure, where surface states forming inside the bulk gap correspond to one- 
dimensional bands indexed by total angular momentum. When a half-integer magnetic flux pierces 
the nanowire, we find a strongly correlated helical Luttinger liquid topologically protected against 
weak disorder. We describe how transport experiments can detect this state. 

PACS numbers: 71.10.Pm, 73.23.-b, 73.63.-b 
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The rich and fascinating physics found in certain spin- 
orbit coupled materials exhibiting the "strong topologi- 
cal insulator" (TI) phase currently attracts an enormous 
amount of attention [l| . In a TI the bulk has a finite gap 
A;, but topologically protected surface modes exist inside 
the gap. Using Bi2Se3, which presently is the reference 
material due to its rather large gap, A& ps 0.3 eV, surface 
probe experiments (ARPES, STM) have provided clear 
evidence for the theoretically predicted masslcss Dirac 
fermion surface state with spin and momentum locked 
together (2j. However, probing the surface state in trans- 
port experiments still poses a major challenge because 
residual bulk charge carriers - either related to disorder 
or due to unintentional intrinsic doping - tend to mask 
the surface contribution even in the cleanest samples so 
far available [||. The surface contribution is easier to 
extract experimentally in thin-film geometries Q or in 
TI nanowires 0, |(| , where the surface-to- volume ratio is 
more advantageous. In the latter case, introduction of 
a magnetic flux $ piercing the nanowire has allowed to 
successfully identify the Aharonov-Bohm effect caused by 
the surface state. 

These recent developments clearly demonstrate the 
need for a comprehensive effective low-energy theory of 
the electronic properties of TI nanowires, which we for- 
mulate here. Very recent work 0, Q has addressed the 
effect of strong disorder for the noninteracting problem. 
We instead consider the weak disorder limit but take 
into account electron-electron (e-e) interactions in a non- 
perturbative way. We obtain the band structure of a 
cylindrical TI nanowire from three different approaches: 

(i) using the low-energy approach of Zhang et al. Q, 

(ii) from the distorted diamond lattice model with spin- 
orbit couplings introduced by Fu et al . II 01) . and (iii) 
using a surface Dirac fermion theory [11] . Taken to- 
gether, these calculations draw a consistent picture for 
the surface states inside the bulk gap, even for very thin 
nanowires: a one-dimensional (ID) electron waveguide 
with modes indexed by the half-integer total angular mo- 
mentum j is formed, where each mode contains a right- 
and a left-mover. The spin direction is always tangential 
to the surface and perpendicular to the momentum. For 



integer flux $ (in units of the flux quantum), we have 
an even number of massive ID Dirac fermion modes, un- 
like the case of carbon nanotubes [ll, ■ This allows 
for impurity backscattering, and with e-e interactions 
one has standard disordered Luttinger liquid (LL) behav- 
ior [IH[lij| , where the SU(2) spin symmetry is broken. 
The case of half-integer $ is more intriguing. Here an 
emergent time reversal symmetry (TRS) for the surface 
states results in an odd number of modes topologically 
protected against weak disorder. With interactions this 
yields a helical Luttinger liquid. In the simplest single- 
mode case, the spin polarization of a right (left) mover 
has a counter-clockwise (clockwise) orientation around 
the waist of the cylinder. The helical LL has been de- 
scribed previously [TH, [I?} as edge mode of the 2D "quan- 
tum spin Hall" (QSH) topological insulator realized in 
HgTe/CdTc quantum well structures [l8|. However, it 
has been difficult to reveal the QSH helical LL experimen- 
tally, since usually the edges living on opposite bound- 
aries both contribute. While more complicated setups 
involving junctions of different QSH systems have been 
suggested [lji[, the situation is unique for a TI nanowire 
at half-integer $: the fermion doubling theorem [l|| is 
circumvented and, effectively, just one QSH edge can be 
realized. This simpler realization of a helical LL should 
allow for clear signatures in transport experiments. 

Let us start with the band structure of a clean nonin- 
teracting cylindrical nanowire for $ = 0. First, we em- 
ploy the low-energy approach of Zhang et al. 0] where, 
expanding up to order fc 2 in momentum around a suit- 
able symmetry point, e.g., the T point in Bi2Se3, the bulk 
TI Hamiltonian consistent with TRS plus inversion and 
rotation symmetry has the form 
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The z direction defines the anisotropy axis, eo{k) 
C + Dik 2 z + D 2 k\, M(k) = M + B x k 2 z + B 2 k\, and 



h 2 — h 2 



k 2 We use Pauli matrices er for spin and 



t for parity (orbital) space; cto and To denote the iden- 
tity. The TI phase is realized for M$B\ : 2 < 0, and we 
take parameters for Bi2Se3 as quoted in Ref. @. For 



FIG. 1: (Color online) Band structure of a TI nanowire with 
R = 15 nm obtained by numerical diagonalization of Eq. 
Points refer to bulk states, lines to surface states. Inset: 
Density (p) (dashed red) and spin density [(s^,): blue, (s z ): 
black curve] vs radial coordinate for the right-moving state 
(fc,j) = (0.02 A" 1 , 1/2). 



a nanowire along the e z direction [we use cylindrical 
coordinates with unit vectors e r = (cos 4>, sin <p, 0) and 
= (— sin 0, cos <f), 0), and put h = 1], rotation sym- 
metry in the xy plane implies conservation of total an- 
gular momentum J z = —id$ + <J z /2, with half-integer 
eigenvalues j. For nanowire radius R, we require the 
wavefunction to vanish at the boundary r = R, which is 
automatically ensured by__ expanding in the orthonormal 
set of radial functions 



,(r < R) 
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where 7 m „ is the nth zero of the Bessel function J rn 
with integer m. For given (k = k z ,j), we express Hz 
in the basis \n,a,r), where a = ± (r = ±) denotes the 
eigenvalue of a z (r z ) and the associated radial function 
is u_j_ CT /2,n( r )- Some algebra gives 

H z \nar) = (e (k) + M{k)r)\nuT) + Arfaln, a, -t>(2) 

2lA 2 \ - 7j+cr/2,n'7j-o-/2,n , , 
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with the substitution fcj_ — > lj~a/2,nl R nl e o(^) an d 
M (k). Numerical diagonalization is then straightforward 
and yields topologically protected surface modes. A typi- 
cal band structure and spin (particle) density profiles are 
shown in Fig. [T] Evidently all surface modes have a fi- 
nite gap. States with (k,j) and {—k, — j) form a Kramers 
degenerate pair, and for given k, the ±j states are degen- 
erate but have opposite s z spin polarization. We observe 
that the expectation values of the spin density operators 
s^^^ = je^ ri2 • cr only depend on the radial coordinate 
r. Since then (s r ) = — ^(s^) = 0, spin is always oriented 
tangential to the surface. Moreover, we find that the spin 
direction always encloses the angle r\ = it/2 with the mo- 
mentum k = ke z + (j/R)e<f,. For large \k\, a right (left) 



FIG. 2: (Color online) Same as Fig. [T] but using the tight- 
binding model on a diamond lattice, see Eq. ©. We take e z 
along the (111) axis, R — 5a, A so = t, and St = 0.28t. The 
size of the unit cell along this direction is d ce ii = y/3a. Inset: 
as in Fig. [T]but for fcdceii = 2.9. 



moving surface state then has counter-clockwise (clock- 
wise) spin polarization (s^) > ((s</>) < 0). 

More microscopically, a TI nanowire can be described 
by a tight-binding model for the electronic states in a 
diamond lattice with spin-orbit coupling A so (Toj . 



r- E 4(-[4><4])^ 

((id)) 



(3) 

where a is the lattice spacing. To reach the TI phase, a 
distortion — > tij + St in the nearest-neighbor hopping 
is introduced along the (111) direction. The A so term in- 
volves second neighbors and depends on the two nearest- 
neighbor vectors d 1,2 connecting them. After choosing 
an axis direction (e z ), the wire is formed by all lattice 
sites located within radius R. A typical band structure 
for a nanowire with R = 5a and e z in the (111) direc- 
tion is shown in Fig. [5J The surface states are again 
characterized by a finite gap, and spin or particle densi- 
ties are qualitatively consistent with those shown in the 
inset of Fig. [TJ The radius dependence of the lowest sur- 
face state gap, A S (R). obtained under both approaches is 
compared in Fig. [31 where we set the parameter — 2t + St 
in H t b, which is half the gap in the (111) direction, equal 
to Mq in Eq. ([T]). Agreement between both models at 
large R is reached by adjusting a to 2.8 nm. We see 
that even for very thin nanowires, the analytical predic- 
tion A s = i>2 /R, see Eq. below, agrees very well with 
the tight-binding result, while the low-energy model @ 
gives deviations when R < 5 nm. It is worth mention- 
ing that even though a k ■ p expansion of H t b around 
the L = -(1, 1, 1) point does not match completely with 
Eq. ((T|), the main features of the surface states are equiv- 
alent in both descriptions. We have checked that similar 
results are obtained when e z points along other crystal- 
lographic directions. 

Both the dispersion relation and the spin texture of the 
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FIG. 3: (Color online) Numerical results for the lowest surface 
state gap A s vs nanowire radius R, obtained from the low- 
energy approach @ [black circles] and from the tight-binding 
model ([3| [red squares]. The analytical prediction v%/R, see 
Eq. (0, is given as blue dashed curve. Inset: $ dependence 
of A s for R = 5a from the tight-binding approach, where the 
flux is introduced via Peierls phases. 



dispersion ([5]) comes from a Berry phase n due to spin- 
surface locking [ll|. While scattering between Kramers 
pairs, (k,j + $) <R- (—k, — (j + $)), is forbidden since the 
states (JSJ) have zero overlap, backscattering (fc — > —k) for 
fixed j is allowed, i.e., potential scattering (disorder) is 
relevant. A non- integer flux $ lifts the degeneracy, and 
for half-integer <£>, the mass appearing in Eq. ([5]) van- 
ishes for the special band j = — <f>. This feature also ap- 
pears in the tight-binding calculation, see inset in Fig. [3] 
Spin-conserving single-particle backscattering processes 
are then forbidden, and weak disorder has no effect (2l| . 
When the chemical potential p is inside the bulk gap, we 
are thus guaranteed to have an odd number of modes. 

In the remainder, we focus on half-integer 3>. For sim- 
plicity, we consider p < A s and sufficiently weak inter- 
actions, where only the single mode j = — <fr needs to be 
retained in a low-energy effective theory. Moreover, we 
assume p > so that Umklapp c-e scattering can also be 
neglected, cf. Ref. [l(|. Using the spinors (jH), the surface 
electron operator ^(z, </>) is expanded in terms of slowly 
varying chiral ID fermions ip r =±{z). 



surface states found under these two approaches are well 
reproduced by a model of 2D massless Dirac fermions 
wrapped onto the cylinder surface, under the condition 
that the spin is tangential to the surface and perpendic- 
ular to the momentum (77 = 7r/2), cf. Ref. [xTj [. To match 
the above numerical results, we also need to take into ac- 
count anisotropy, since there are different Fermi velocities 
Vi_ t 2 along the e z ^ directions. By supplementing Eq. ^ 
with boundary conditions describing a flat 2D surface in 
the xz plane, we find v h o = A 1;2 ^l - (D2/B2) 2 . Taking 
the parameters of Ref. [!| for Bi2Se3, V2/V1 ~ 2. With 
i] = 7r/2, the surface Hamiltonian takes the form [Tl| 



1 



r=± 



1p r (z) 



r 



(7) 



with the Fermi momentum fc_F = pjv\. The stan- 
dard bosonization approach [bfj expresses ipr{z) — 
(27r£) -1 / 2 exp[iv / 7r((^ + r6)] in terms of conjugate phase 
fields f(z) and 0(z), where the surface layer width £ 
is the short distance cutoff for the ID continuum de- 
scription. The noninteracting Hamiltonian is Hq = 
^ / dz[(d z ip) 2 + (d z 9) 2 ]. The density operator, p(z, 4>) = 
v&tv]/, is then equal to the ID density, d z 9 / ' ^pn. Similarly, 
the spin density operators and s z are reduced to a ID 
form, 



surf 



-ia z <p/2 



v\ka y 



V2 , 



,i<7 z c/>/2 



(4) 

where we added the dimensionless flux parameter $. We 
note that 4> may include not only the orbital magnetic 
field, but also a Zeeman field or an exchange-coupled 
magnetization due to a nearby magnet (for fields along 
e z ). The dispersion relation implied by Eq. (|4|) is 



E k , j>± = ±\v 2 k 2 + 



R 2 



(5) 



where the ± sign refers to conduction and valence band, 
respectively. The corresponding eigenstate is 



U k,J,± 



±ijl 



l k,j,± 



U k ,j,± 



y/2E k , ji± (E k j,± + (j + W)V2/R) 



(6) 



For integer $, all bands are doubly degenerate and have 
a gap > A s = V2/R- The mass term in the relativistic 



s z 



1 



V(z, 



J{z) 

■^cos[2k F z + 2^9{z)} 



(8) 



We observe that equals the ID current density, J{z) = 
d z (p/^/TT 1 reflecting spin- momentum locking. There are 
no 2k f oscillatory terms in p nor in s^. On the other 
hand, no "slow" terms exist for s z , and we always have 
(s z ) = 0. 

We now include e-e interactions, assuming that no 
metallic gates are nearby. Similar to the nanotube case, 
apart from a hard-core part, the main contribution to the 
(surface-projected) potential can be modeled by [Tj| 



U{r 



+ £ 2 + 4i? 2 sin 2 [(</>- 4>')/2] 



where k takes into account the dielectric constant of 
the substrate and of the insulating interior part of the 
nanowire. Inserting the field operator ([7|) into the gen- 
eral second-quantized interaction Hamiltonian yields the 
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ID expression 

He ~ e = hj dzdz ' v{ - z ~ zl) d ' e (^w (9) 

with the effective ID potential V(z) = 
(27r)- 1 d(f> U(z,(f>). The explicit form of V is 
given in Ref. [ 1 3| and has the Fourier transform 
V(q) ~ (2e 2 /K)[0.51 - \n\qR\] for \q\R < 1. Hard- 
core interaction terms give an additional contribution 
b J dz[(d z 9) 2 — (d z (f>) 2 ] to the Hamiltonian, where 
b depends on microscopic details. Since b stays 
marginal under renormalization group transformations, 
the logarithmic singularity in V(q), caused by the 
long-ranged Coulomb tail, is expected to dominate 
in practice. Approximating q « 2n/L for nanowire 
length L, we obtain the single-mode helical LL fl6| . 
ffhix = jjdz [K(d z ^) 2 + K-\d z 6f] , where v = v x /K 
and 

K = 1 (10) 

V'l + ^ r (ln[L/(2^)] + 0.51) 

with K = 1 without interactions. It is st raig htforward to 
generalize these expressions to b ^ (T3.ll9j| . For L/Rm 
1000, we find the strongly correlated value K w 0.4 to 
0.5 from Eq. We note that K « 0.53 to 0.9 for 

the QSH edge in HgTe/CdTe wells. 

The helical LL state in a TI nanowire can be identified 
through several experimentally observable signatures. 
First, we note that the equal-time spatial correlations of p 
and = J decay as \z\~ 2 . While (s z ) ~Q remains valid 
for arbitrary K, s z correlations show a slow algebraic 
power-law decay, (s z (z)s z (0)) oc cos(2k pz)\z\~ 2K . For 



K < 1, we therefore find an ordering tendency towards 
spin density wave (SDW) formation, where spins are ori- 
ented along the nanowire axis e z . Within the standard 
classification of ID systems [l4|, the helical LL in a TI 
nanowire is thus in a SDW phase. As a consequence, the 
Ruderman-Kittel interaction among magnetic impurities 
mediated by such a nanowire is extremely anisotropic (see 
also Ref. [22| ) and decays only with a slow power law. At 
the same time, the absence of Ikp terms in the density 
operator implies that no charge density wave (CDW) cor- 
relations develop at all. Furthermore, the superconduct- 
ing order parameter describing singlet Cooper pairing is 
0(z,4>) ~ e % ^ip+(z)ip_ which implies a fast power- 
law decay cx \z\~ 2 / K . The angular e" 1 ^ dependence comes 
from the spin structure in Eq. ([7]) and causes an addi- 
tional strong suppression of the proximity effect when 
bulk superconductors are in contact to the nanowire. For 
normal-state metallic electrodes, in a two-terminal geom- 
etry, the conductance is G — e 2 /h (independent of K) 
when the contacts are ideal. However, non-ideal contacts 
cause a typical temperature-dependent decrease of G(T) 
at low temperatures due to the well-known power-law 
suppression of the tunneling density of states [19| . More- 
over, in contrast to a spin-polarized LL (which also has 
G = e 2 /h for ideal contacts), spin plays an essential role 
here. This could be easily seen in the presence of mag- 
netic impurities. In particular, the Kondo effect can take 
place, where theoretical predictions for G(T) exist [23| 
and directly apply. To conclude, we are confident that 
the helical LL will soon allow for its clear experimental 
identification in topological insulator nanowircs. 
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